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I . Introduction 

In studies of the propagation of finite amplitude sound 

I 

waves in liquids, recourse is generally made to the so-called non- 
j linear parameter D/A, which can be defined by an expansion of 
pressure changes in terms of density changes in the liquid for an 
isentropic process. If we use a Taylor expansion for this relation 



(?£) <P-pJ + (~ §) 



(P-Po>‘ 



° '°P'o,s'' 'O' 'ap 2 o,s 2 



( 1 ) 



where p,p Q , are the instantaneous and hydrostatic pressures, P, P Q , 
the instantaneous and equilibrium densities, where the derivatives 



are taken at P = P Q and constant entropy, and limit ourselves to 
terms through quadratic, we can rewrite the equation as 



P = P c + A 



P-P, 



P-ft ^ 2 



a - Ic-^) 



( 2 ) 



Since the first derivatives is the square of the infinitesimal 

2 2 

sound velocity, c Q , A = p Q c , and we can write the ratio 



B _ JL(JDl) 
A " c 2 a P 2 



(3) 



o,s 



1 ) 

2 ) 

3) 



There are three principal method of determining B/A; 
Finite Amplitude Determination, 

Static Determination, and 
Thermodynamic Determination. 



Finite Amplitude Determination 

This method makes use of the harmonic distortion in a 
propagated wave. The amplitude of the second harmonic is compared 
to the amplitude of the first harmonic. There are several 
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different experimental approaches to the problem. 

i 

Shklovskaya-Kordi^ used an acoustic filter to filter out the 
second harmonic and a quartz receiver to obtain data. Adler and 
Hiedemann^ used the same arrangement except that the measurements 
were taken by determining the variation of the index of refraction 
of the liquid instead of using a receiver. Various other experi- 
mental techniques are discussed by L.K. Zarembo and V.A. Krasil'iUkov/' 
Some problems associated with the above method are that 
the filter employed may interact with the system, and if it passes 
second harmonic waves then it passes the fourth harmonic also. 
Reflections may also be a problem especially if an acoustic 
receiver is used. Theoretically, the equations concerned with 
finite amplitude are difficult to manipulate and, in some cases, 
cannot be solved without approximations on the non-linear effects. 

Static Determination 



A static determination of B/A requires precise measure- 
ments of density as a function of pressure at constant temperature. 
The data is analyzed for the slope and derivative of the slope of 
specific volume with respect to pressure. This information will 

i 

provide the non-linear parameter of the isothermal equation of state , 
The difference between this parameter and the isentropic non-linear 

' s 

parameter has been investigated by Coppens. While the difference 
is quite small, accurate measurements of the thermodynamic 
parameters needed to calculate this difference are almost non- 

i- 

i existent. 

.} 

The main difficulty with this method of determining B/A 

3 

is that very accurate data are required as the two derivatives are 
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l 

about equal and are subtracted from each other. 



t Thermodynamic Method 

This method involves measuring the velocity of sound 
land its derivatives with respect to pressure and temperature . 

i 

l Because of the simplicity of this measurement compared to the other 
determinations , the thermodynamic method was chosen for use in this 

!• 

research. 



There are two devices which measure the velocity of sound 
with the accuracy desired, the velocimeter and the interferometer. 
The velocimeter was chosen, as the interferometer, because of its 
movable parts, would require elaborate seals to contain the 
pressures required. 

To establish the thermodynamic method, we rewrite 
equation (3). 



f = ^<-2% = = 2p^C&> 



c 0 o,s 9(5 °j s ° 0 o,s 



the differential dS can be written: 



dS = (||) p dT + (f ) T dP 

Using Maxwell's relations, we recast the above equation in the form 

TdS = CpdT - T(||) dP 

For an isentropic process 
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Now 



d0 = ( t ) T,S dP + ( ^>P,S dT + < "^ ) T , p dS 

= ( -of4,S dP + % ( f ) P,S ( -i>P,S dP + ( f > T,P dS 

(_c>c.) _ (^C\ + PVT / Son 

op's op ^ T , S C p ^6T'p,S 



where 



8= i(~) 

p V V vxP,S 



Substituting the above expression in equation (3)> we obtain 

- = 2p c [(— ) + (J.G) i 

A o o L ^ep ; T,S C p ' dT ; P,S J 



= 2p C (~) 



where 



00 op T.s 



= (B/A)' + (B/A)" 



C p = pV * 



( CO) 

— ' pn' O Q 

C ? Ox * ,s 



(3h) 



Thus we have adequation where the components can be experimentally 
determined. 

’ 

The assumption that sound transmission is isentropic has 
been investigated by Gol'dberg.^ He showed that if the absorption 
of sound is small over a wavelength (i.e., a/k « 1) and 



WK ( 1 — — ) 

dp 

72 



« 1 



men 



po v c 



(0p/0s) o (s-s o ) 



c, r 

(1 -r 1 ) 



(Gp/Cp) (P-P ) Dc .2 

1 S O r<-w<w 
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where a is the absorption coefficient, k the wave number, u/2tc the 
frequency, ;; the coefficient of thermal conductivity, c v (c_ ) ) tho 
heat capacity for constant volume (pressure), P the density, and 
c the velocity of sound. In most liquids the absorption coefficient 
a is sufficiently small to satisfy Gol'dberg's requirements and in 
the present case experimental observation indicates that a is 
sufficiently small here also. 

) Cz v 

«K(1 - c-> * 

5 — 2- ~ 2 x 10~ 5 « 1 

P° v c 

for the liquids investigated. Coppens^ also shows that higher 

I 

order terms and the cross terms of the complete Taylor expansion are 
negligible. 

I 

I 



! 
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II. Experimental Equipment and Procedure 

A . Method of measuring sound velocity . 

The method used for the measurement of sound velocity in 
this research was the velocimeter method of Greenspan and Tschiegg^ 
as adapted by Coppens.? The sample was contained in a cylindrical 
tube each end of which was terminated by a 5 Mc/s quartz. < 
crystal for the transmission and reception. of sound pulses. A 
knowledge of the repetition rate of the sound pulses is sufficient 
to determine the time required for a sound pulse to travel the 
length of the velocimeter. When the length of the velocimeter is 
accurately known the velocity can be determined. 

! 

The principle of the timing method can be understood by 
means of Fig. 1. 



Initial 

Pulse 




transmitting receiving 

crystal (a) t=0 crystal 




transmitting receiving 

crystal (b) t=T crystal 



Second 

Pulse 





(d) t=3T 



Fig 
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Sound pulses are introduced at t = 0 and t = 2T and are detected by 

4 

I 

a second transducer at the opposite end of the velocimeter. The 
tine T can be adjusted so that the two pulses will arrive at the 

1 

receiving transducer at the same time (Fig. 1). The repetition 

|, 

rate was adjusted so that the two pulses add for a maximum 
amplitude. 

Looking at Fig. 1, a and c, we see that the repetition 
rate is the reciprocal of the time required for the pulse to travel 
twice the length of the velocimeter, so that c = ~ * 

The equipment required to measure the velocity of sound 
as a function of temperature and pressure divides naturally into 
four systems. 

The pressure system provides and regulates pressures up 

E o 2000 psi which can accurately be measured to within 2 psi. 

The temperature control and measurement system provides 
temperature range from 4-6°C to 523°C for the sample when the 
ressure vessel is placed inside a high temperature oven. The 
easurement system consists of a platinum resistance thermometer 
connected to a Mueller bridge. A spotlight galvanometer is used to 
balance the bridge. 

The electronic system produces the electrical pulses that 
are changed into the sound pulses by a transducer, accurately 
measures the repetition frequency of the pulses, and displays the 
received pulses on an oscilloscope to allow a correct superposition 
v0 be made. 



•1 
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The pressure vessel and components consist of the 
velocimeter and associated equipment, and a container to allow 
pressure to be put on the sample. 

| i 

I 

1. Pressure system : 

A pressure line from the vessel leads to a Heise-Bourdon 
tube gauge and to a dry nitrogen tank which provides pressures up 
to 2200 psi. The pressure in the vessel is adjusted by a standard 
I high pressure regulator. The system also has a provision for 
1 allowing air evacuation of the pressure vessel by a vacuum pump. 

2. Temperature control and measurement : 

The platinum resistance thermometer is placed in the 
receptacle provided in the base of the pressure vessel. This 
position inside the metal tends to average the temperature 

fluctuations of the oven and also permits a temperature reading 

5 

near the sample. The galvanometer used with the Mueller bridge 
)is a high sensitivity spotlight type. It has an internal resistance 
of 78 ohms, sensitivity of 0.00 5// amp /mm, and a period of *+.3 sec. 

3. Electronic system : 

A sine wave audio-oscillator provides the signal which 
eventually (when appropriately modified) excites the quartz trans- 
ducer. The frequency (repetition rate) of the oscillator is 
accurately counted by an electronic counter. The sine wave passes 
through a band pass filter into special electronic pulse-shaping 
j units. The band pass filter reduces electrical line interference 
‘ hat affects the <Hs; • .y o*> v .10 oscii ioscope. The power supply for 
^-he special electronic pulse-shaping units was a constant voltage 
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power supply. Direct current output was regulated to better than 
; 1 /j over a range of 200 to 325 volts. An isolated AC voltage of 
6.3 volts is also available for tube heater circuits. The sine 
| wave is first led to a Schmidt trigger where a train of rectangular 
pulses is produced at the audio oscillator frequency. Tho shape 

of the various signals are shown in Fig. 2. The signal then goes 

J 

to blocking oscillator #1 which produces spikes at the audio 
oscillator frequency. These spikes go to a blocking oscillator 
counter. The input signal is then used to excite an RC network 
which does not allow the device to pass any other input signals 
for a certain time interval. The tine interval can be changed by 
external adjustments. This signal is then led to a monostable 
multivibrator. The output here consists of a series of rectangular 
pulses. The frequency of the audio oscillator is a multiple of 
this pulse frequency. This signal and the signal from blocking 
oscillator #1 are now led to a gate circuit. The signal from the 
multivibrator acts as a switch to allow a certain number of spikes 
(usually 2) from blocking oscillator #1 to go to the final blocking 
oscillator. The interval between these two spikes is 2T and is 
determined by the frequency of the audio oscillator. The final 
blocking oscillator shapes the signal again, and the signal is sent 
to the transmitting transducer in the velocimeter. The signal from 
'the receiving transducer is then amplified in the wide band chain 
■amplifier and sent to the oscilloscope for display. The trigger 
input for the oscilloscope is taken from blocking oscillator #2. 



1 

a 










1 
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b . 2l ”"-ss(?l and components : 

The main, component of the pressure vessel is the veloci- 
vor . Ti.e quartz c ,y line in the velocimeter is manufactured 

I 

fro:,, dmersill II quartz (an optical quality fused quartz)[see 
Fig. 3~ •> It is made in the form of a hollow cylinder with an 
outer diameter of about 1 3/16 in. and an inner diameter of 5/8 in. 
There is a hole •£- in. in diameter perpendicular to the drilled 
portion which is used for filling and emptying the delay line. 

The ends are polished while the other surfaces have a ground finish. 

1 The maximum and minimum lengths are diametrically opposed and 
measure 1*998995 and 1.998930 in,, respectively. The maximum 
deviations of the ends from perfect planes are 30 and 50 x 10“^ in. 
The length of the delay line is greatest toward its center (around 
the bore). The linear thermal expansion coefficient a for quartz 
is 7*2 x 10"*' 7 /°C~10/2 while the isothermal bulk compressibility K T 

O 

is 18.6 x 10“°/psi+7$* These data were supplied by the manufacturer.; 

The transducers are 1 3/18 in. X-cut quartz discs with a 
3/8 in. concentric disc of gold over chromium plated on one side. 

The metallic sample provides the electrode on the other side. All 
surfaces are polished and the fundamental frequency is designed to 
be approximately 5 Kc/sec. ' 

The pressure vessel itself and the metallic components 
were machined from austinitic annealed stainless steel (type 321, 
3V7K This type does not react with the metals used. 

A pressure inlet is provided in the side of the vessel 



with a line leading to the pressure generating apparatus, 



An 



•opening for the pi. 



resistance capsule is provided in the 



Sleeve 




Cap with 
Transducer 
and Gasket 



Delay 
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' solid material in the bottom of the pressure vessel. 

The quartz delay line, cushioned by gaskets, is placed 
^n a metallic sleeve (see Fig. 3)* Caps containing electrical 
connections hold the transducers tight] y against the delay line 
with the aid of a gasket to obtain clamped boundary conditions, 
i cup with a hole in its center is positioned over the opening in 
the side of the delay line so that the metal can run into the 
cavity upon melting. Since bismuth expands upon freezing, the 
above equipment must be able to rotate so that the metal can be 
poured from the openings in the delay line into a receptacle under- 
jneath the above parts. This assembly is attached to the pressure 
vessel's cover by electrical leads from the cover to the caps. 

An annealed copper gasket is used to seal the pressure vessel. 

The cover is then bolted to the pressure vessel. 

i. 

| 

5. List of envtnmer.t : 

j 

Pressure vessel and components 
I Composition gaskets 

Mfr: Armstrong Cork Company 

Asbestos-nitrile rubber composition An-890 Accopac 
Pressure gauge 

I Heise-Bourdon Tube gauge temperature compensated; 
rated to 2000 psi 
High pressure regulator 

Ufr: Linde Division of Union Carbide Corporation 

High pressure fittings and piping 

Vajc 3T * J ^xCcln ^ u cnp Cl 



I 

c 



r 




__O3052 

__0584 

J03841 

05842 




XAA 

CM. 



0586 

p586l 

0587 

0388 

05881 
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High temperature oven 

Mfr: American Instrument Company 



'•'•Inn T- T • r-> 


Tu. pera 


ture Bench Mode Oven: 


Cu w*v_» xCjj it O t 


' r — 3 ^ 1 1 0 




Temperature 


range : 


46 *C to 523*0 



Maximum temperature gradient: 2°C 

Temperature regulation: 0. 3*C 
Platinum resistance thermometer 

Mfr: Rosenount Engineering Company, Model 162D 

Mueller resistance bridge 

Mfr: Gray Instrument Company, Model S-1003 

High sensitivity spotlight galvanometer 

Mfr: Rubicon Division of Minneapolis-Honeywell 

Regulator Company, Model 3432 
Wide band chain amplifier 

Mfr: Spencer, Kennedy Laboratories, Incorporated, 

Model #202D, Bandwidth 100 kc to 200 Mc/sec. 
with a voltage gain of 20 db + 1 db. 
Oscilloscope 

Mfr: Fairchild Dumont, Type 766-H/P2B 

X - Plug in unit 

Type 74-13 delaying sweep with precision delay 
X - plug in unit 

Type 7 6-02 dual trace 

Frequency response: DC to 25 mc/sec. 

Amp', lav a e : to 5 mv/cm. 

<r Mlo os.; ill r 

Mfr: IIcvlco<; Packard Company, Model: 200 CD 
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Electronic counter 

Mfr: K'.'jct'j Pac’.card Company, Model 523 DR 

Variable band pass filter 

Mfr: Xrohn-Eitc In strum on.; Company, Model 320-A 

Constant voltage power supply 

Mfr: Associated Specialties Company, Model 3 

Quartz delay line 

Mfr: Anersil Quartz Division of Englehard Industries 

Incorporated 
Transducers 

Mfr: Valpey Corporation 



C . s'.-** -'.nr- emeri mental apparatus 

me i * ± a m used ras a sample that had been used 

'previously in measurements of 3/A, but since it had been used in 

ji 

; the presence of a nitrogen atmosphere in a quartz and stainless 
steel container, it was considered reusable. The original metal 
had been obtained from the Indium Corporation of America with a 
guaranteed purity of 99.999/®* The bismuth was obtained from 
Fisher Chemical Company (Cat. Mo. B-319) with an estimated purity 
of 99*98$. Analysis of the metal as certified by the supplier 
jlisted the following metallic impurities: 



Arsenic 


(As) 


0.0000^ 


Copper 


(Cu) 


0.002 % 


Iron 


(Fe) 


0*001 $ 


- D 


(?b) 


0.01 % 


..r/cir*ony 


' (Do) 


r, c, c * 

J# UU 
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Before any data for the alloys were taken, sound velocity 
.measurements were obtained for indium which agreed within 0.6 m/sec. 
i 0 . 00 ^ n/ooe-°C to that of Coupons . 

Velocity measurements wore made at four concentrations 
| of bismuth. According to the 3i-In phase diagram, there are two 
eutectic and one peritectic points (at 34$, 66$, and 48$ Bi 
respectively) „° Data were recorded at these points and at a 17$ Bi 
concentration. The latter point was studied in order to check on 
the approach of the alloy's characteristics to that of indium. 

Some measurements were also taken for bismuth. 

Before any measurements \. r ere made, the experimental 
apparatus was first thoroughly cleaned with dilute nitric acid, 

I water, benzene, and acetone rinses. 

The first experimental sample was prepared by cleaning 
' the metal as above except that the acid rinse was with a much 
weaker solution because the surface of the metal tended to oxidize. 

I 

The metals were then -weighed in order to determine and obtain the 
desired percentage alloy. They were then put into a beaker and 
heated under paraffin oil to about 300°C. The mixture was held at 
this point for about five minutes and was stirred vigorously. The 
alloy was then allowed to cool while the stirring continued. The 
freezing point of the alloy was determined roughly and checked with 
the phase diagram. The alloy was then cleaned wi th the above 
mentioned rinses and placed in the cup above the delay line. 

Since indium is a fairly high priced metal, a technique 
f successive di.1ution was used. The samples were tested in order 
of increasing a; • r-.tw of bismuth. After each sample had been 

tad. nn npprnpri a t-o ar.. mt of it was cleaned and weighed for the 



It 
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next sam. c. The amount of bismuth necessary for the next 
concentration was put into a beater with the required amount of 
he previous saapje and the see unique described above was again 
followed. The final sample was analyze^ by Jarrell-Ash Company 
and found to have a concentration of 66. 33$ Bi. The computed 
percent bismuth was 67.00$. Therefore, the maximum error in any 
i of the concentration values is considered to be about 1$. The 
computed concentrations were 17.00$ Bi (hereafter designated as 
17$ 3i. ) , 33 .84$ Bi (hereafter designated as 34$ Bi.), >+7.66$ Bi 
(hereafter designated as 48$ Bi.), and 67.00$ Bi (hereafter 
designated as 66$ Bi.). 

In order to get the acoustical signal into and out of the 
(Sample, a wetting agent was used. A silicone, Dow-Coming 200 
fluid (cs 350, 12,500) served this purpose quite satisfactorily. 

The transducers were thoroughly cleaned with acid, benzene, and 
acetone rinses. They were then coated with the wetting agent. 

The coating had to be very thin in order not to influence the 
velocity data. In all cases the coating was thin enough so that 

interference patterns (Newton's rings) were seen. The components 

l 

/ere then assembled and the whole assembly v/as placed in the oven, 
lectrical, pressure, and thermometric connections were made. Then 
the oven v/as brought slowly (to prevent warping of the pressure 
vessel) to the desired temperature. Air v/as evacuated from the 
vessel and nitrogen admitted. This v/as done to insure that the 



i. 



■ogen atmosphere and therefore its surface 



-“1 



13 



D . " 1 i by i 2 ' rl " m c- 1 1 

After the sample had melted (as indicated hy the 
appearance of a transmitted signal on the oscilloscope), several 
hours were a 2 lowed to p: ss before any readings were taken. Since 

i, 

ch~ thermometer was not in she sample, equilibrium concitions 
would only bo determined by checking the change in repetition rate. 

A minimum of 3 hours ras allowed to elapse between pressure readings 
and S hours sc tween temperature readings. It was found that after 
1 hour, the repetition rate had reached 90 % of its final value, 
after 2 hours 99 % of its final value; so that in a three hour 
period, the eapuilibrium error was well within r.m.s. reading error, 
-ha c.'ght hour period wa s allowed to elapse for a temperature 
change. Data taken for a combination temperature and pressure 
change indicated that in a four hour period the repetition rate 
was 99 % of the final rate. 

2. Measurement of repetition rate and tenmerature 

The readings were taken by adjusting the repetition 
rate so that the superposition of the two pulses at t = 3T (see 
fig. 1) was a maximum. This procedure insured that the time 
between the pulses was the time required for one pulse to travel 
to one end and return. The wave packet at t = T was also displayed 
on the oscilloscope to aid in determining correct superposition. 

The frequency of the audio oscillator was then determined from the 
electronic counter. Temperature measurements were taken at this 
1 time also. Groups of ten frequencies at ' to 7 temperature 
:• endings v w when together uu til he average of two groups of 



- , eqae..cy 



within 0.1 cycle. The more consistent 
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roan-ngs were cnen usee lor c. 



Results of using the above methods were satisfactory for 



•27 



t 



- * o vh 



apoe 



O ~ \ f . 

o J CuuoCu vuG 



made to obtain 
r six days abov 



ed. 


, . Hoi 


/ever , 


or 


only 


three 



bismuth. The high molting point of bismuth 
"wetting” agent to decompose. Four attempts 
data. The last attempt was the most successful, 
c molting temperature, a signal of correct form 
the amplitude kept decreasing each day so that 
different temperatures were obtained. 



Ii 
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1 . S-^np^norn t.~ r~" c t 

The tran: eel on cocff'I dent for the transducer in 
velocimeior is near aero f •_ r th air backed sice and near 
unity for the s ar.pl o sice of the crystal. These values indicate 
energy is transferred very easily from the transducer and to the 
sc rule and bade. A very short ■’.rave packet is therefore produced 
as she energy given the transducer by the electrical pulse is 
jsickly dissipated into the sample. At the receiving transducer, 
interference may arise as some components of the wave packet arrive 
• t the front of the transducer at the same time when some other 
components have reflected from the rear of the transducer. The 
interference effect can be computed by a straightforward analysis 
when the input is assumed to be a delta function.- 1 ^ The result of 
this calculation is shown for indium ^.n Fig. 4 ^ An actual 
: waveform for a sample composed of 66 % Bi is shown superimposed on 



the computed values for pure indium. The curves almost coincide 

6 

because the acoustic impedance for bismuth is 1,67 x 10 while the 
acoustic impedance for indium is 1.62 x 10^ cgs units. Looking at 
the experimental curves we see the difference in amplitude to be 

j 

slight between the first and second pulse which shovs that one of 
Gol’dberg's required conditions (a/.: <<1) is satisfied. Notice 
that the first i-g- cycles of the calculated and experimental pulse 
at t = 3? are completely missing, while the next half cycle is very 
smell* This effect is the result of the interference mentioned 
above. Figure s- can also be h ought of as the picture of the 

■ppw ” " :fcr . ;-g the second pulse intro- 

>nin a eerrea t ropusiaion r..,te, the pulses should be 



^.4- 'i r> j* v 

6 o — \s. a - 
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A. Delay line correction 

Since the length of the delay line changes with 

temperature and pressure, an error is introduced thereby in the 

determination of the velocity. It can be shown that the length 

3 

change is governed by the equation 



l = £ o [l+cc(t-20°C)][l - ~ K T p] 



OO 



III. Corrections and Ca li brations 

The equipment's design provides a high degree of accuracy 
in measuring relative sound velocity in the delay line. However, 

I errors in the absolute determination of the velocity of sound may 

I arise from two different sources. The instruments such as the 
thermometer, pressure gauge, and electronic counter should be 

I 

| calibrated with the absolute values of their respective quantities 

I 

so that absolute valurp of the sound velocity can be measured as a 
function of absolute quantities. The fact that measurements are 
taking place in a delay line instead of free space also influences 
these measurements. While some of these conditions are systematic 
and will not greatly affect the accuracy of the derivatives with 
respect to pressure and temperature, an attempt is made to correct 
for all possible errors. 



The maximum change in length of the delay line amounted to 0.0017 cm 
Which caused a change of 0.5 m/sec in velocity. Because of the 
uncertainty in the •• .ues of a and Krji , the maximum uncertainty due 



to this correction was about 0.003 m/sec. 
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superimposed as shown. The only pulses that could be added are the 
first half of the third cycle. The only error introduced by this 
superposition is due to waveguide effects. However, since these 
pulses are small, the changes in the sura of the amplitude is also 
small. Better accuracy can be obtained in finding a maximum 
i amplitude by superimposing pulses of larger magnitude. Therefore, 
the pulses with the arrows drawn above were superimposed. A 
correction to the velocity of sound must now be made since the 
•ocond pulse enters the delay line a certain time later in order 
to superimpose at the new position. Since the delay is one wave- 
length over a distance of 2-t, the correction to the sound velocity 
is 

( 5 ) 



Ac = fx 



Ac varied from 5.225 to 9. 664 m/sec, with an accountable uncertainty 
of 0.003 m/sec. One sound velocity measurement was taken using the 
correct superposition and compared to one using the above super- 
position. Differences in velocity amounted to 9.454 m/sec, all of 
which could be accounted for by equation ( 5 ) except for 0.210 m/sec. 
Experimental uncertainties amounted to 0.177 m/sec, so the unaccount- 
able uncertainty this "incorrect" superposition introduces amounts 
to about 0.04 m/sec. This uncertainty is probably due to waveguide 
effects . 

C . V. r avegulde effect correction 

Waveguide effects account for the next largest 
correction. When sound travels along a waveguide, the various 



•equency components travel at different speeds, depending on the 



24 



frequency. This produces a dispersive effect. The leading edge of 
i each cycle does travel at the free field velocity, but since 

measurements were taken at maximum amplitude, this effect will 

I 

cause a systematic error. The amount of this error can be 

i 

calculated for a particular arrangement.^ The error in the 
measured velocity of sound due to waveguide effects varied between 
0.160 and O.38O m/sec. The uncertainty introduced by this correc- 

J 

tion is negligible. 



D. Viscosity correction 

Viscous forces in a medium tend to slow the passage 
of sound waves. Kinsler and Frey state the expression for the 
velocity change Ac as^-3 



Ac = c [- 



:] 






( 6 ) 



ifttpva 2 

where p/p is kinematic viscosity, v is frequency of sound, and a is 
the radius of waveguide. This correction ranges from 0.100 to 
0.049 m/sec. Values of kinematic viscosity were found only for 
indium^I^Bi^, and bismuth^ Values required for other concentra- 
tions were calculated from above values. Assigning an error of 50 % 
to the value of kinematic viscosity used gives an estimated 
juncertainty of 0.004 m/sec here. Another possible source of error 
in the interpolated data may be the variation of viscosity with 
pressure,^ but since no data could be found on its magnitude in 

liquid metals (in organic liquids, the viscosity roughly doubles 

| p 

•for a pressure change from 1 to 2000 kg/cm ), no estimate was made. 
;Since the viscosity correction is small and a 50 % error was assumed 

I 

jin the value of kinematic viscosity, any variation due to pressure 



i r~ ~ ~ * 
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should lie within these limits. 

I I 

F. . Tempe r ature errors 

The maximum temperature error is dependent on the 
accuracy of the calibration of the thermometer because the fluctua- 
tions in the reading of the bridge were within the accuracy of the 
calibration of the platinum thermometer. 

F. Calibration 

In order to determine the absolute velocity of sound 

ll 

as a function of temperature and pressure, the accuracy of the 

i 

thermometer, pressure gauge, and electronic counter was ohecked 

j 

against absolute values. 

The platinum resistance thermometer was calibrated by 
Rosemount Engineering Company and found to be accurate to within 
+0.015°C from -183°C to +350°C. 

The Mueller Bridge has a guaranteed accuracy of ,+0.005$ 
+0.00005 ohms at a temperature of 23 ±1°C. 

The Heise-Bourdon Gauge was calibrated with a dead weight 
piston gauge and found to be accurate within +2 psi. 

The electronic counter was calibrated with WWV and 
adjusted for correct readings. It has a stability of two counts 
in one million per week and is accurate to +1 count. 




i 

1 
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IV. Data 

A • £LBJ2 pi e calculatio n f or velocity of sound . 

A calculation is shown in detail below for one of the 

I' 

points of the Z ] '% Bi sample. The experimentally determined data 

i 

< 

are : 

Repetition rate 20.76668 +0.00077 kc. 

Temperature 77. 48 +0.02°C 

Pressure 0 +. 2 psi 

I 

Actually the temperature is not determined directly, but is 

! 

obtained from a table where resistance values are listed opposite 
the temperatures. 



Corrections for changes in the length of delay line: 



2 1= 2-t o [l+a(t-20 o C) ]pL - j ] 

= (10.15^752 cm) (1.0000414) (1) 
= 10.155172 cm. 



(4) 



Hence 



c = 2 U = (10.155172 cm) (20.76668) kc. 
= 2108.892 m/sec. 



Correction for superposition: 




= (2108.892) 2 

(10.155172) ( 5*io 6 ) 



= 8.759 m/sec. 



1 1 



'i 
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Correction for waveguide effects: 



1 = IQ'MKzJ. 
w 2 TT^ 



a 2 /X 



x - Distance pulse travels 
a - Radius of delay line 
X - wavelength in medium 



(p/w) - xc — -r- = X4.35*10" 6 )c = 0.922171 

(0.6>+)(5*10 b ) 



© = cot 



J _ — 1 • 



10003 radians 



1 + 



-g"(P/w) 



2 

, J(P/w) + 9 



wt 0 = ; 7 = 1*195793 

1 + ^(p/w) -3 



(p/w) ^ = 0.307390 

x=5 

WtQ = 1A28036 



© = l.H -1709 radians 



wt 



c = 



3600 

Correction for viscosity: 

It 

h — { 

O ' 1 

xp 



S - 



= 0.331 m/sec. 



Ac = c f = c [ 1 - *" - ^ - — ] = 0.0V2 m/sec. ( 6 ) 

u+xpva^ 1 +x(5*10 6 )(0.6 l +) 



Velocity of sound corrected for changes in length of delay line 

2108.892 m/sec 



1 . 

ij 

i, Correction Cor superposition 

|j 

Correction. Cor wave;*,.: I de effect 



\ on 



V > j I - ■» 



8.759 m/sec 
O. 33 I m/sec 
0.042 m/sec 
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Error due to uncertainties in temperature 

and pressure corrections of delay line: 0.003 m/sec. 



Error due to uncertainties in superposition 

correction: 0.003 m/sec. 



Error due to \mcertainties in viscosity 



correction: 



0.00*+ in/sec. 



Error due to root mean square deviations in 



Corrected sound velocity 



obtaining repetition rates: 



0.078 m/sec. 
2117.362 +0.088 m/sec. 



The experimentally determined sound velocities of the 



117$ Bi concentrations are given in Table I. The results of the 
alculations for each set of concentrations with the same pressures 
were fitted by a least squares method to polynomial equations by 
the Brown University IBM 7070 computer. The results of the 
computation for the sound velocity of the 17$ Bi sample are also 
given in Table I below the experimental values of c. Polynomial 
curves of degree one through three were calculated. The linear 
least square fit resulted in a smaller root mean square deviation 
from the experimental points of the 17$ Bi alloy than the second 
md third degree polynomial equations while the second degree 
polynomial least square fit resulted in smaller r.m.s. deviations 

from the experimental points of the 3*+$ Bi, *+8$ Bi, and 66$ Bi alloy 

/ 

than the first and third degree polynomials. However, there was 
still quite a large curve fitting error, even larger than experi- 

f 

mental error in some cases. 

Table II, which tabulates sound velocity as a function 
of temperature and pressure was computed from the polynomial 
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equations in the case of the alloys. For bismuth, the points 
Civon are the experimental results. 



slues for the velocity of sound in a Y?% Bi-In liquid alloy corrected for superposition 
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i 8. Velocity of Sound Tables 

TABLE II 

Voikos of the velocity of sound (me tors/sec . ) as determined from 
the polynomial fitted to the experimental points. 



17% Bi (by weight) -In alloy 







Pressxire (Psi) 






To mo 

°c 


0 


500 


1000 


1500 


2000 


! 125 


2222.419 


2223. 066 


2223.595 


2224.012 


2224.518 


jl'-o 


2217-525 


2218.140 


2218.663 


2219.113 


2219-608 


161 


2212.630 


2213-214 


2213.732 


2214.214 


2214.698 


ll80 


2207-464 


2208.015 


2208.527 


2209.043 


2209.515 


202 


2201.482 


2201.994 


2202.500 


2203.055 


2203.513 


l . - 


34% Bi (by we 


ight)-In alloy 




77 


2117-485 


2117.972 


2118.558 


2119.045 


2119.532 


98 


2113.580 


2114.118 


2114.641 


2115.151 


2115.665 


125 


2108.281 


2108.869 


2109.364 


2109.882 


2110.404 


143 


2104.573 


2105.186 


2105.696 


2106.204 


2106.717 


.161 


2100.727 


2101.357 


2101.907 


2102.394 


2102.885 
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48% Bi (by weight) -In alloy 




98 


2013.542 


2014.024 


2014.584 


2015.176 


2015.711 


125 


2010.789 


2011.219 


2011.835 


2012.431 


2012.933 


143 


2008.535 


2008.953 


2009.607 


2010.182 


2010.676 


161 


2005.945 


2006. 370 


2007. 061 


2007.599 


2008.094 


180 


2002.847 


2003.300 


2004.030 


2004.510 


2005.017 


194 


2000.325 


2000.812 


2001.356 


2001.995 


2002.517 


66% Bi (by weight) -In alloy 
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1896.378 
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1893.524 
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1895.398 
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! 180 


1886.174 


1886.787 


1887.463 
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i!l9*f 
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1884.414 


1885.180 


1885.745 


1886.347 
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1880.517 


1881.029 


1881.633 
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Pressure (Psi) 



rn 

ip 

0 
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500 


1000 


1500 


2000 


Bismuth 


313 


1639. >+33 


1640.174 


1640.943 


1641.866 


1642.542 


323 


1637.768 


1638.552 


1639. 206 


1639.927 


1640.612 


341 


1637.05-0 


1637.647 


1638.191 


1638.991 




357 


1633-933 












Since the 


errors for 


all but the 


deviations 


when obtain- 


ing repetition rates 


5 are very small, they 


can be considered the 



sane for all points. The r.m.s. error varied from about 
0,030 m/sec to a maximum of 0,075 m/sec. The curve fitting error 
also varied widely from a negligible error up to 0,137 m/sec with 
an average of 0,52 m/sec. 

From earlier discussion: 

Error due to uncertainties in 

temperature and pressure corrections: 0,003 m/sec. 

Error due to uncertainties in super- 
position corrections: 0,003 m/sec. 

Error due to uncertainties in viscosity 
correction: 0.004 m/sec. 

Error due to root mean square deviations 
in obtaining repetition rates: 

Average 0.075 m/sec. Maximum 0.116 m/sec. 

Total error of experimental points 






Average 0.085 m/sec. l.iaximum 



Maximum curve fitting error 



Average curve fitting error 



0.126 m/sec. 
0.137 m/sec. 
0.052 m/sec. 
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C . Compu t ation of B/A 

After the equations of sound velocity at constant 

» 

pressure were determined as a function of temperature, equations 
for sound velocity at constant temperature were determined as a 
function of pressure. These curves were plotted and found to have 
a linear dependence. Graph V shows the curve with the maximum 
slope ana C?raph VI the curve with minimum slope. 

Beside the value of c and its derivatives, the computation 
of 3/A required three experimentally determined values. These are 
density, thermal expansion coefficient, and heat capacity of the 
samples . 



Values for the density of liquid In-Bi mixtures could not 

oe found. A density curve vs. $Bi(by weight) at room temperature 

1 8 I 

(i.e* for the solid) showed an approximately linear relation. 

4 1 7 

^Equations to find density as a function of temperature ' were 

i 

‘extrapolated and density values computed using a linear relation. 
*Kow r ever , since bismuth expands upon freezing and the freezing point 
of ala. the alloys tests was below that of bismuth, an attempt was 
made to determine the density and change of density as a function 

of temperature for the 66$Bi sample. The data obtained agreed 

i 

roughly with the calculated values. 

The thermal expansion coefficient was not found for the 
alloys, but since indium and bismuth each have the same value for^ 
the coefficient and the term containing this quantity is by far the 
j smaller of the two making up B/A (Eq.3a), negligible error should be 
introduced here. 

Heat capacity (C,) values for the alloys also were 
jhnavai lable. A linear dependence variation between the values 
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T found for bismuth and indium was therefore used .17 Any error in this 
value should have a very small effect as it is also in the smaller 
term of B/A. 

The results of the calculations for density and heat 



capacity are listed in Table IV. Values of B/A were computed from 
the values listed in Table IV according to equation (3), and are 
tabulated in the table. 
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The error in determining B/A for the alloys is mainly the 
uncertainty of the values of (B/A) 1 since (B/A)" has an order of 
magnitude equal to the uncertainty of (B/A) 1 . Due to the estimated 
error in the density values (+0.5 gm) , an uncertainty of 0.15 is 
calculated for B/A. The uncertainties in velocity are very small 

las the max:' .urn uncertainty introduced to B/A is 0.00008. The 
uncertainty in the slope of c vs. pressure affects the accuracy of 
B/A greatly since the slope of this line is very small. Errors in 
the data points on the order of +0.1 m/sec. will cause an 
uncertainty of 0.3 in B/A. The dependence of c vs. p appears to be 
linear well within experimental error, so this assumption was used 
in the computing uncertainties above. The use of the average value 

I 

of r.m.s. errors and curve fitting errors resulted in an average 

j 

uncertainty to B/A of 0.218. The use of the maximum values 
resulted in ah uncertainty of 0.58. 

I! 



The calculation for the uncertainty for bismuth follows 

the same procedure as above, with a much smaller uncertainty due 

. 

to density and a much larger uncertainty due to the slope of c vs. 
p as data at only three temperature points were used. 



Errors for alloys 



1 Error 

I 


in 


(3/A) 


O 


to uncertainties 


of isity 


0.15 




1 Error 

I 


in 


(3/A) 


due 


to uncertainties 


of c 


0.0008 




!• 

Error 


in 


(3/A) 


due 


uncertainties of 


(3c/5p) .j.(av) 


0.218 


(max. ) 

0.^33 


Total 


uncertainty 


in (3/A) 


(av) 


0.37 


(max. ) 
0.58 
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Errors for Bismuth 

Srror in (B/A) due to uncertainties of density 

• Error in (3/A) due to uncertainties of c 

Error in (3/A) due to uncertainties of (3c/3p). 

\i ' 

1 

i Total uncertainty in (B/A) for bismuth 



0.09 

0.0004 

0.96 

1.2 
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V. Cone - ; us ions 



As far as could be 

Sj 

work done in determining B/A 

{ 

measurements for pure indium 

as a function of temperature 

I. 

results with the appropriate 



determined there had been no previous 
for indium-bismuth alloys. B/A 
and velocity measurements for bismuth 
have been reported. The available 
results of this work are given below. 



Temp ( C) 

/* 

( 6c/Cp) xiO”°(cm/dyne-sec) 



(3c/2t) 

P 

c (m/sec) 
.(B/A) ! 
(B/A)” 
.(B/A) 

Temp (°C) 
271 
285 
289 
318 



(cm/sec- C) 



Indium 

(Coppens) 

160 

1.47 

-28.5+0.6 
2313±0.6 
4,78 
-0.23 
4 • 5+5$ 



17 % Bi-In Alloy 
(This work) 

161 

1.49 

- 27.4 to. 0 . 

2212 to. a 
5.02 
-0.26 
4.8 



Bismuth 

Velocity (m/sec) 

1635+15 

1663+5 

1666 

1639+2 



(3c/9t) (m/sec°C) (Source) 

-0.5 (20) 

( 22 ) 

0 ( 21 ) 

-O./fc (This work) 



The results of this work appear to approach the values 
:-r indium as the concentration of indium is increased (see also 
lections IV-3 and D) . The velocity of bismuth found in this 

‘research agrees with Kloppa 1 s value if (Oc/Qt)^ found in this work 

! ' 

us*;-’ to ektrapoiate to the tempera tur o used. Otherwise the 
■ •' : 1 v,--> re lies h • : ween the two r.-.nges of reported values. 
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It is possible that the sound velocity measurements 
reported hero were made under non-equilibrium conditions or condi- 

I 

tions which were not considered when mal&g applicable corrections. 
Tor about six days after the sample had melted, the signal on the 
oscilloscope indicated poor acoustical contact of the sample with 
the transducers. Then the signal changed to one indicating good 



acoustical contact. However, the voltage amplitude of the signal 
started to decrease slowly, so that after five more days an 
accurate superposition could no longer be made. Upon disassembly 
of the pressure vessel and components, the transducers were found 
to be cracked. Bismuth had leaked out of the delay line and was 
found in the sleeve between the transducer and cap at one end of 
the velocimeter. The reason for loss of amplitude of the signal 
was probably due to this leakage of bismuth. The bismuth collect- 
ing behind the transducers may have affected some of the earlier 
assumptions used in computing corrections such as the assumptions 
of air-backed transducers, a filled delay line, and a fundamental 
resonating frequency of 5 Mc/s. The bismuth probably leaked 
through a crack in one of the transducers caused by thermal 



[stresses set up by the high temperatures. Due to these uncertainties 
n error in velocity measurements of 2 m/s is assigned. A larger 
error is not considered reasonable because the waveform observed 
jjdid agree very closely with the calculated waveform (see Section 
ilIII-3), but since the set of velocity measurements was not completely 
consistent, the above estimate of error was made. 

Hudriavtsev derived an expression for the velocity of 

• pi 

ction of concentration of the components. 



;ounc in mixtures a 



CL ^ Cv 



I 

•fl 
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Assuming an ideal mixture, we have 



U = ,uU-i + (l-;i)U, 



(7) 



where U is the internal energy and ;i is the molar fraction of the 
first component. Also 

v = + (l-;a)v 2 

. where v is the volume. Differentiating equation ( 7 ) twice with 
respect to volume at constant entropy and noting that 



•l 



(I&l) = 

v „ 2 s 

S v v 



I where M is the molecular weight and c is sound velocity, 

1 

Kudriavtsev obtained 

o 2 hp 2 

c = y ” c-, + ( 1 -a) zr 31 

M M d 

where M is the average molecular weight. A plot of the above 

values of sound velocity is shown on Graph 7 along with the experi- 

24 

mental results of Coppens and this work. The maximum difference 
)be tween values of c is 

The reason the c vs. temperature dependence for three of 
| the alloys has a varying slope is not known but may be due to the 
difference in slope for this dependence in the pure metals or may 

| 

vindicate a change in the slope of temperature dependence of bismuth 

i 

’beyond the range investigated. 

The 3/A data for the indium -bismuth system are shown on 
Graph 8 , The values increase regularly with concentration with the 

C 

^exception of the 4vv Bi alloy. The compound In 0 Bi is formed, and 
a peri tootle was found by IVrotti in this area of concentration. 



concentration. 
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No other unusual characteristics have been reported. 

The non-linear parameter of the equation of state 

P = p(?,s) bas not been derived from first principles, nor has 

' 

even a rough guide been established for predicting values. Thus 

} Si 

far, experimental values have been obtained for some metallic 
elements, alloys, water, and various organic fluids. Reported 

♦I 

> S 

values of B/A range from 2.7 to about 10.5. 

There are a few physical constants that appear to be 
related to B/A. The metals provide a range of B/A values from 2.7 
to 7*8. Physical constants that follow this trend are 1/c, 

1/latent heat of vaporization, and density (see Graphs). However, 

i when water and organic liquids were considered, the only physical 

constant to follow the trend was 1/c. The fact that 1/c appears 

to be proportional to B/A is somewhat surprising as in equation 3a, 

2 

it is directly proportional to B/A. Acoustic impedance (pc), pc , 

j |j 

and B do not follow this trend. Table V lists some of these 
parameters. 

The effect of temperature on B/A may be important in 
determining what physical mechanisms determine the degree of non- 
linearity. Water changes by a factor of 1.9 over the temperature 
range 0°-100°c , while ethanol's maximum variation in B/A is 0.6. 

The fact that there are less parameters that follow the 

I trend of B/A for organic liquids is not considered unusual, because 
the organic molecules are much more complex so that combinations of 
other effects may be more important. 




TABLE V 51 

Physical Constant?; of Some Substances 



J 

Material 

! 

1 


Temp 

(°cj 


B/A 


Velocity Latent 

(m/seo ) heat f 

vaporiza- 

tion 

( col/gm? 


Density 

(gm/cm 8 ) 


Isother- 

mal 

Compres- 

sibility 

(cm^/dym'/j 


Sodium 


110 


2.7 b 


2521 b 


1005* 


0.93 a 


21 . 0 h 


S Tin 

\ 


240 


4.4 b 


2470 b 


57 3 a 


6.95 a 


3.1 h 


i 

Indium 


160 


4.5 b 


2313 b 


468 a 


7.03 a 


3*2 h ! 


Water 


10 


4.6 C 


1447° 


590 g 


1 . 00 s 


5.0 s 


Bismuth 


318 


7.1 


1639 


204 a 


10 . 03 a 


1 

1 


Mercury- 


30 


7.8 b 


l440 b 


70 a 


13-52 a 


0.4 i 


Benzene 

1 


30 


9.0 d 


127^ 


94@80° S 


0 . 88 s 


10 . 2 s 


5 

Acetone 


20 


9.2 d 


ll 88 d 


125@56° s 


0.79 s 


12 . 8 s 


i 

Methanol 


0 


9.^ 


ll89 f 


263@65° g 


0.79 s 


10.7 s 


Methy- 

Acetate 


30 


9.7 d 


1128 d 


ll4@58° s 


0.92 d 


— 


Diethy- 

lamine 


30 


10.3 d 


109^ 


91@58° g 


0.71 s 




n-Propanol 

j 


0 


10 . 5 * 


129 5 f 


164@97° s 


0 . 80 s 


10 . 5 s 


Cthanol 


0 


10. 5 d 


1232 f 


20 l >@78 g 


0.79 s 


9.9 s 



a. 


Ref. 17 


f. 


Ref. 28 


b. 


Ref. 19 


g- 


Ref. 29 


c. 


Ref. 2 5 


h. 


Ref. 20 


d. 


Ref. 26 


i. 


Ref. 13 


e. 


Ref. 27 










-Propanoi . Ethanol 
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GRAPH 9 
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